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Abstract. We study a class of S"^ Gowdy vacuum models with a regular past Cauchy 
horizon which we call smooth Gowdy symmetric generalized Taub-NUT solutions. 
In particular, we prove existence of such solutions by formulating a singular initial 
value problem with asymptotic data on the past Cauchy horizon. The result of our 
fT^ ■ investigations is that a future Cauchy horizon exists for generic asymptotic data. 

Moreover, we derive an explicit expression for the metric on the future Cauchy horizon 
in terms of the asymptotic data on the past horizon. This complements earlier results 
O , about X §^ Gowdy models. 



PACS numbers: 98.80.Jk, 04.20. Cv, 04.20.Dw 
1. Introduction 



Studies of cosmological solutions of Einstein's field equations have a long tradition and 
led to astonishing results about our own universe. In particular, observations indicate 
O I that there was a big bang in the distant past, and indeed, the simplest cosmological 
models, namely the Friedmann solutions for reasonable matter fields, predict precisely 
this behavior. The question arises whether such curvature singularities occur for generic 
■ solutions of Einstein's field equations or if the strong symmetry assumptions underlying 
^ I the Friedmann models are necessary for this. The Hawking-Penrose singularity theorems 
[T7] shed some light on this question. They predict incompleteness of causal geodesies 
in a wide class of situations. However, the information about the reason for the 
incompleteness provided by these theorems is very limited, and it is indeed not true 
that it is always caused by a geometric singularity. 

Let us restrict all of our investigations to vacuum with a vanishing cosmological 
constant and to four spacetime dimensions. The corresponding Cauchy problem for 
Einstein's field equations is well-posed and leads to the notion of the maximal globally 
hyperbolic development (MGHD) of a given Cauchy data set [U [35]. The example of 
the Taub solution [37] shows that incompleteness of causal geodesies, as predicted by 
the singularity theorems, signals a different kind of phenomenon where it is possible 
to extend the MGHD. The extended solutions are not globally hyperbolic. There 
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exist closed causal curves and indeed, there are several non-equivalent extensions. 
This unexpected property has caused an ongoing debate in the literature. Are such 
pathological phenomena a typical feature of Einstein's theory of gravity - in which case 
it could not be considered as a "proper" physical theory - or do such phenomena only 
occur under very strong and special assumptions, for example the high symmetry of the 
Taub solutions? An interesting and intensively debated hypothesis in this context is 
the strong cosmic censorship conjecture whose widely accepted formulation was given 
for the first time in [11], based on ideas by Eardly and Moncrief |2^ and Penrose |32j . 
More details and references can be found in IM] . This conjecture states that the 
MGHD of generic Cauchy data is inextendibleo- Roughly speaking, this implies that, in 
the generic situation, incompleteness of causal geodesies is indeed caused by a geometric 
singularity in some sense. 

In his effort to generalize the family of Taub solutions and hence to show that there 
is a large (but still non-generic) class of solutions of the field equations with similar 
"undesired" properties, Moncrief defines the family of generalized Taub-NUT spacetimes 
in He is able to prove an existence result under an analyticity assumption. Our 
motivation in this paper is two- fold. First, we want to extend the existence result to 
the smooth cas^ and formulate it as a singular initial value problem with "asymptotic 
data" on the Cauchy horizon. However, in this paper, we restrict ourselves to the 
case with Gowdy symmetry. Secondly, we want to study the global behavior of such 
solutions. By means of so-called soliton methods, it turns out that in the generic case, 
the existence of a past Cauchy horizon implies the existence of a future Cauchy horizon, 
at least under certain topological assumptions on the horizons discussed later. The 
reader should compare these results to the x §^ case in [T9] . 

The paper is organized as follows. In Section |2] we discuss some background 
material, in particular the symmetry reduction introduced by Geroch which is necessary 
to be able to write the metric and the field equations in a useful way. Section |3]is devoted 
to the definition and discussion of our class of smooth Gowdy symmetric generalized 
Taub-NUT solutions based on Moncrief's earlier mentioned class. The existence and 
uniqueness theory and the corresponding singular initial value problem is considered in 
Section 13.31 The basic ingredients for these investigations are Fuchsian methods which 
we describe briefly in the appendix. The next main part is Section H] where we discuss the 
global-in-time properties of smooth Gowdy symmetric generalized Taub-NUT solutions. 
We finish the paper with conclusions in Section |5l 

^ This conjecture, of course, only makes sense if one is able to give a precise and reasonable meaning 
to the terms "generic" and "inextendible" . At this stage, however, this has not been found for general 
situations. See e.g. [55] . 

We use the term "smooth" for infinitely difFerentiable objects, as opposed to "analytic" for which 
the Taylor series converges in addition. 
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2. Geometric preliminaries 

2.1. Symmetry reduction by Geroch 

We briefly present here the symmetry reduction introduced by Geroch in [T3]. Let 
M = M X if be an oriented and time-oriented globally hyperbolic 4-dimensional 
spacetime endowed with a metric g^^ of signature (— , +, +, +), a global time-function 
t and a Cauchy surface H. We denote the volume form of (/^^ by e^^pa and the 
hypersurfaces given by t = to for any constant to by Ht^. Each Ht^ is homeomorphic to 
H. 

Now, let ^ be a smooth spacelik^ Killing vector field which is tangent to the 
hypersurfaces Hf and set 

The twist 1-form of ^ is 

and V is the covariant derivative compatible with g. The field ^ is hypersurface forming 
if and only if fi^ = 0, which will, however, not be assumed in the following. We define 
the "3-metric" 

on M, and, by raising indices with the inverse of g, we define also h'^i, and h'^'^ on M. 
The first of these tensors is the projector to the space orthogonal to ^ in TpM, p G M. 
From the volume form e^j^p^r of (7, we furthermore introduce the tensor 

1 

Let be any 1-form. One defines the derivative operator D as 



'^pup ■ r—'^pupa<, ■ 



Note that at this stage we are only interested in local patches of M. Then, the 
flow generated by ^ induces a map vr from M to the space of orbits S, i.e. vr maps 
every p G M to the (locally) uniquely determined integral curve of C, starting at p. The 
requirement that vr is a smooth map induces a differentiable structure on S, and hence 
S can be considered as a smooth manifold. Using C^h = and h{^,-) = 0, Geroch 
shows that there is a unique metric on S which pulls back to /i^j/ along vr. We call this 
metric on S again /i^,^. The same can be done for h^^^, and h^'^, which are henceforth 
considered as objects on the quotient space S. The first is just the identity operator, 
while the second is the inverse of the 3-metric h^^, on S. We can proceed in the same 
way with the function A, the 1-form the tensor e^^i^p and the covariant derivative 
operator D^, which can henceforth be considered as objects on S. Then ef^^p becomes 
the volume form of /i^^^ and the covariant derivative operator compatible with /i^,^. 

^ Note that Geroch considers the case of a timehke Kilhng vector field. 
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Geroch shows that Einstein's vacuum field equations on (M, g) imply that the 1-form 
fl is closed, dQ = 0. Our local considerations allow us to introduce the twist potential 
u so that Q = du. The basic quantities A, u and h^i, on S completely characterize the 
geometry of {M,g). 

Now Geroch introduces a conformal rescaling 

We refer to the associated covariant derivative operator as D, Ricci tensor as S^j^i, etc. 
He shows that the vacuum field equations for (M, g) (and certain geometric identities) 
are equivalent to the following set of equations 

b'x = J (^b^xb.x - b^cob.co^ , (i) 
b^co = ^b^xb^co, (2) 

A 

s^u =^{b^xbA + b^ujb,ujy (3) 

These equations are the Euler-Lagrange equations of the Lagrangian density 



C= \ -deth 



5 + ^ {b^'xb^x + b^'ujb^uj 



Hence, the equations can be interpreted as 3-dimensional gravity on S coupled to a 
wave map u : S ^ H where "H is the 2-dimensional hyperbolic space represented by the 
components (A,^). 

We point out that the quotient manifold 5* can be identified with the surfaces 
orthogonal to ^ if and only if ^ is hypersurface forming, i.e. fi^ = 0. However, in general 
this is not the case and the manifold S cannot be interpreted as a submanifold of M. 

2.2. Spacetimes of spatial 3-sphere topology with isometry group U{1) or U{1) x f/(l) 

We specialize the previous general discussion to the case M = M x S'^. Now, the 
Cauchy surfaces are H = E>^. We think of as the submanifold of determined by 
x1 + X2 + xl + x1 = 1. The Euler coordinates of §^ are written as 

xi = COS - COS Ai, X2 = COS - sm Ai, 

X3 = sm - cos A2, X4 = sm - sm A2, 

where 9 G (0, tt) and Ai, A2 G (0, 2tt). Clearly, these coordinates break down at the poles 
9 = and vr. We will also use the coordinates {9,pi,p2) determined by 

Ai=:(pi + P2)/2, A2=: (pi-p2)/2. (4) 

Note that the coordinate fields and dp^ are smooth non-vanishing vector fields on 
S^, while the fields dx-^ and dx^ vanish at certain places. Note that these fields can be 
characterized geometrically (without making reference to coordinates) in terms of left- 
and right-invariant vector fields with respect to the standard action of SU{2) on see 
for example [H |5] . 
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We specialize now to the case of a metric g with a spacehke Kilhng vector field ^ 
which generates a smooth effective action of the group U{1) on M such that ^ is tangent 
to the level sets of the time function t. In particular, the integral curves of ^ are closed 
in Ht for each t. Moreover, we assume that H, and hence Ht for each t, is a Hopf bundle 
generated by ^. This means that the map vr, defined in the previous section, is a Hopf 
map from M = M x §^ to = M x S^. In order to give an explicit representation of this 
map, we can introduce, on Hf for every t, either Euler coordinates or the coordinates 
Xi, X2, X3, X4 of the embedding into R^. There is no loss of generality in assuming 
^ = dp-^ because we can apply a diffeomorphism from to itself to achieve this for 
each t. However, from the point of view of the initial value problem, the assumption 
C, = dp^ is not compatible with all coordinate gauges; see below. In any case, this allows 
to represent the map vr as 

TT : M X §^ ^ M X §^ 

(t,xi,X2,X3,X4) {t,yi,y2,yz) 
= (t, 2{xiX3 + X2X4), 2{x2X3 - xiX4),xl + xj - xl - xl) (5) 
= (t, sin 9 cos p2 , sin 9 sin p2 , cos 9) . 

Here we consider §^ as the submanifold of determined by + + ?/| = 1. When 
we introduce standard polar coordinates on S^, namely 

?/i = sin cos 0, ?/2 = sin 1? sin 0, 2/3 = cos "(9, (6) 

then we have that vr is simply 

{t,9,p,,p2) ^ {t,^,<j)) = {t,9,p2). 

In particular, the push-forward of dp^ to M x §^ along tt vanishes and the push-forward 
of dp2 equals the coordinate vector field on 

We point out that Eq. yields a global representation of the map vr. Since both S 
and M are simply connected, the twist potential oj is defined globally. The vector fields 
dp-^ and dp^ will play a major role in the following discussion. The fields dx-^ and dx^, on 
the other hand, cannot be used for Geroch's reduction directly because they vanish at 
either ^ = or ^ = vr. 

Let us now consider a global smooth effective action of the grou F@?7(l)xf/(1) with 
a second spatial Killing vector field rj commuting with ^, under the same assumptions 
as before. We caiifl suppose ^ = dp-^ and rj = dp^, since it can be shown that all smooth 
effective actions of U{1) x U{1) on S'^ are equivalent in the sense that any other smooth 
effective action equals the previous action after possibly applying a diffeomorphism of 
into itself and an automorphism of U{1) x U{1) to itself [lO]. This group action 
degenerates, in the sense that the group orbits become 1- dimensional, precisely at ^ = 
and at ^ = TT. When ^ = 0, the field 8x2 vanishes, while dxj^ vanishes at ^ = vr. Note 

^ This is the direct product of the group U{1) with itself. 

But this assumption leads to further restrictions on the possible choices of coordinate gauges for the 
initial value problem, see below. 
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that the fields dp^ and dp^ never vanish, but both become parallel at = and 9 = tt. 
The action of the group U{1) discussed in the previous section can be associated with 
the action of the subgroup {e} x U{1) of U{1) x U{1). 

We consider Geroch's symmetry reduction only with respect to ^ and the 
corresponding projection map vr according to Eq. ([5]). Since ^ and rj commute, we 
are allowed, in principle, to apply the symmetry reduction successively one after the 
other for both fields. However, the result is then not a smooth manifold, but rather a 
manifold with boundary. In order to avoid the corresponding complications, we perform 
the reduction only with respect to ^ and obtain the orbit manifold = M x §^ with 
the 3-metric h with Killing field rj. Here, the push-forward of rj along tt to S* is again 
denoted by rj. Recall that the push-forward of C, is 0. Now, according to Geroch [TS] 
and Gowdy [T6], Einstein's vacuum field equations imply that the twist quantities of 
the f/(l) X f/(l)-action on M, 

vanish. The geometrical interpretation is that the 2-space orthogonal to the 2-space 
spanned by ^ and 77 in M is integrable to a 2-surface everywhere. This allows us to 
make the following general ansatz for the metric g on M, 

g = gAEdx^dx"" + R [e\dpi + Qdpa)' + e'^pa'] , (7) 

where A, B = 0,1 label coordinates on the submanifolds orthogonal to the Killing vector 
fields; the metric gAB is so far unspecified. The functions R, L and Q only depend on 
t and 6, i.e. are constant along the Killing vector fields. It also follows that A and u, 
as objects on S, are constant along 77. Moreover, rj is a hypersurface orthogonal vector 
field in 5* (but not necessarily in M). We can compute 

X = Re^, (8) 

from Eq. ([7]), and 

dttu = -i?e2V|det((?cD)|/^5^Q, deoj = -Re^^l detigcD)W^dAQ. (9) 
Finally, the 3-metric is 

R^ 

h = gABdx^dx^ + i^e'^dpa^ = gAsdx^dx'^ + — dpa^ (10) 

A 

Chrusciel [9j defines the notion of smooth generic vacuum Cauchy data on a 
spacelike hypersurface of {M,g). He shows that there exist a smooth functioE0 M, 
a constant Rq > 0, and functions Q and L as above, on (0, vr) x which are constant 
along ^ and 77 such that the spacetime (M, g) with M = (0, vr) x and g of the form 
Eq. ([ID with 

i? = i?osintsin^, (^ab) = e^^diag (-1, 1), (11) 

can be isometrically embedded into the maximal globally hyperbolic development of such 
data. We will assume that Chrusciel's genericity condition is satisfied and henceforth 

The function M is not to be confused with the manifold M. 



Smooth Gowdy symmetric generalized Taub-NUT solutions 



7 



identify {M,g) with [M,g). The question whether {M,g) is isometric to the maximal 
globally hyperbolic development in general is, in some sense, the main content of this 
paper. One calls such a spacetime (M, g) a Gowdy spacetime and the time coordinate 
t is called areal time. We will assume that t & I := (0,7r). 

2. 3. Smoothness conditions for the metric components 

Of particular importance for our further studies is the behavior of the metric coefficients 
of h at the poles of the 2-sphere. We know that h is smooth on / x §^ and invariant 
along Tj. Since 7] is purely spatial, it is hence a rotationally invariant metric. Moreover, vj 
is hypersurface orthogonal with respect to h. We thus consider the standard embedding 
of §^ into in a neighborhood of the north pole 1/2, 1/3) = (0, 0, 1). We introduce a 
local coordinate patch (^1,^/2) there which is related to spherical coordinates (■»?, 0) by 

yi = sin '& cos 0, 7/2 = sin i? sin 0. 
A general smooth metric / has the form 

I = hidyl + 2ii2dyidy2 + h2<^yl- 
We want to study the implications of two conditions now: 

(i) The metric is diagonal with respect to the 0)-coordinates (because 77 = is 
hypersurface orthogonal) . 

(ii) The metric is invariant under rotations around the polar axis (because 77 is a Killing 
vector field). 

This easily implies that there are smooth functions F and G which only depend on cos •& 
such that 

h^{-d^<P) = G{cos'd) - F{cos'd)sm^'dcos{26), 
l22{^,(j)) = G{cos^) + F{cos^) sin2^?cos(20), 
/i2(^9,0) = -F(cos^9) sin2^9sin(20). 

With respect to spherical coordinates, the metric then has the form 

I = cosM [G(cosi?) - F(cosi?) sin^-f?] di?^ + sin^-f? [G(cos^) + F(cosi?) sin^t?] d(/)^ 

This representation is clearly only valid in a neighborhood of the north pole, in particular 
it breaks down at the equator due to the cos^'j^-factor. If we choose the functions F and 
G such that 

Z = e^di92 + sinMe2^d(/.2, (12) 

then M and U must be two smooth functions of cosi? alone, so that the following 
condition holds 

gM _ ^2U ^ M{cos^) sin^i? (13) 
for some smooth function M. 
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According to Eq. ffTOj) . the general metric on 5* = / x is of the form 

h = e^'{-dt^ + de^) + sin^e e^^dps', (14) 
in areal gauge given by Eq. ( ITTj) . Hence 



[/ = In i?o + In sin t - i In A = (In /2o - ^ - In sin 6^ + In sin t) /2. (15) 

This is finite everywhere on / x under the previous assumptions because the quantity 
A must be finite and bounded away from zero at the poles for all t G (0, vr), and hence 

In particular, the smoothness condition Eq. f|T3|) translates to 

d2 

e^^ = ^ + M(cos t?) sin^t?. (16) 

A sin 6 

A consequence is that e^"'^ is bounded and non- vanishing at the poles for all t G (0, vr). 

Our choice of {C,,ri} = {dp-^,dp^} as the Gowdy Killing basis has further important 
consequences due to the fact that dp-^ = dp^ at 9 = and dp-^ = —dp^ at 9 = n. At the 
poles we must have g{^,^) = giv^v) = ^9{^:V)- Therefore there must exist a smooth 
function Q which only depends on t and cos 9 so that the function Q in Eq. ([7]) satisfies 

g(t,^) = cos ^ + Q(t, cos ^)sin2^. (17) 

In particular it follows from Eq. that for each t G (0, vr) 

-2 = g(t,7r) -g(t,0) = / Qgd9 = - R-^e-'^^dtud9 = - RX'-'^dtU d9. 

Jo Jo Jo 

2.4- Reparametrizations of the Gowdy orbits 

All of our discussions so far are based on the choice {dpj^jdp^} as the Gowdy Killing 
basis on M. Now we study general reparametrizations of the Gowdy Killing orbits in 
M, i.e. arbitrary bases of the Gowdy Killing algebra. For example, the standard basis 
normally used in the literature is {8x^^,8x2}- 

Let (01,02) £ be coordinates on the Killing orbits so that a Gowdy invariant 
metric has the form analogous to Eq. ([7]), i.e. 

g = e^'{-dt^ + d9^) + R [e^(d0i + Qdcjy^f + e-^d02] . 

We are allowed to reparametrize the orbits by means of constants a, b,c,d E M, so that 

ad — be 

and 

01 = a0i + 602, 02 = C01 + (i02. 

The coordinates t and 9 are not changed. In terms of the new coordinates, we want to 
write the metric as 

g = e^{-dt^ + d9^) + R P(d0i + Qd02)^ + e-^d0^ 
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One finds that 



R 

Q 



\ad — bc\R, 

(a + c(5)V + c^e"^ 

lac? — bc\ 
(a + cQ){b + dQ)e^ + cde-^ 



(18) 
(19) 

(20) 



{a + cQye^ + c^e-^ 

A particularly useful transformation is the inversion, i.e. the interchange of the 
Killing basis fields. Then we have a = d = 0,b = c=l, and hence 

R = R, 



Q 



e^Q^ + e-^ 



Let us now consider a metric in the parametrization (pi, P2) of the Killing orbits as 
given by Eq. ([7]). We pick 0i = pi and (^2 = P2- Now let 0i = Ai and 02 = ^2 and hence 



g = e^\-dt^ + de^) + R 



For this we must choose a = 1, 6 = 1, c 
R = 2R, 

I (1 + Q) V + e-^ 



e^(dAi + QdA2)' + e- 
= 1, d = — 1 from Eq. 



^dA^ 



It follows that 



Q 



-(l-Q2)ei + e- 



(l + Q)2e^ + e-^ ■ 
The inverse of this reparametrization is 

^ - 2' 

^ (i + g)V + e-^ 



(21) 
(22) 

(23) 



^ (l + Q)2e^ + e-^ 
From this and the discussion in Section 12. 3^ we can easily derive the behavior of 
the functions R, L, Q at the poles for t E (0, vr) in areal coordinates, 

R =Rsm9, (24) 
6 

=e^cot-, (25) 

Q =(l-cos^)Q, (26) 
R 



e^(l - cos 6*) + e"^(l + cosi 



(27) 



with smooth functions R, L, Q and M which only depend on t and on cos^. 

Note the following interesting general fact about polarized Gowdy spacetimes. A 
Gowdy metric is called polarized if there exists an everywhere orthogonal basis of Gowdy 
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Killing fields. With respect to this basis, the function Q must hence vanish identically. 
Now, Eq. f[T7|) shows that this can never happen for the Killing basis {dp-^,dp^}, but it 
is possible for the basis {d\-^,d\^} according to Eq. fl26|) . Indeed, one can show that Q 
can only vanish identically for a smooth Gowdy symmetric metric on §^ if the Killing 
basis is chosen such that one of the two fields is proportional to c^a^ and the other to 

3. The class of smooth Gowdy symmetric generalized Taub-NUT solutions 
3.1. The Tauh solutions 

The Tauh solutions were discovered by Taub [37] as a family of cosmological solutions 
of the vacuum field equation with spatial §^-topology. They are a two-parameter family 
of spacetimes 

,2 / 4(1 + r2)^ 2 2^. ^ A 

9 = 1 y — ^ ^ ''^^^ ®uji + uj2® UJ2) + Y^^^ ® ^0 , 



with / > 0, m G M and 
Here, 



V{t) := -4r2 -8— r + 4. 

6 



ui = sin pid9 — cos pi sin 9dp2, 
U2 = cos pid9 + sin pi sin 6'dp2, 
U3 = dpi + cos6'dp2, 

are the standard left-invariant one-form^ on The metric is smooth and globally 
hyperbolic for 

r G (r_, r+), r± := -m/l ± a/1 + m'^/P. 

It was demonstrated for the first time in [31] that these solutions can be continued 
analytically through the apparently singular times t±. These extensions were christened 
Taub-NUT solutions. The extended solutions are not globally hyperbolic; in particular 
there exist closed causal curves. Moreover, there are several non-equivalent analytic 
extensions. 

We easily find that the Taub solutions are Gowdy symmetric (but not polarized), 
and we can bring them to the form Eq. ([7]) with flTT]) 

g = e*^(-dt2 + de^) + R [e^(dpi + gdp2)' + e-^dp2'] . 

For arbitrary parameters / > and m G M, the Taub solutions are then given by 

R = 2lVP + m2 sin t sin 6, 

e^' = /2 ^ (m + V/2 + rn? costf, 

^ The term left-invariance here is understood with respect to the standard action of SU(2) on Note 
that SU(2) is a subgroup of the isometry group of the Taub solutions. 
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R 



Q 




3.2. Generalizations of the Tauh-NUT solutions 

Moncrief [23] introduces the family of generalized Tauh-NUT spacetimes as globally 
hyperbolic spacetimes (0, 6] x with a smooth time function t for a sufficiently small 
6 > 0. The level sets of t are Cauchy surfaces diffeomorphic to S'^. Moncrief wants to 
study the situation where global hyperbolicity breaks down at t = 0, in the sense that 
the spacetimes can be extended through t = as non-globally hyperbolic spacetimes so 
that the points corresponding to t = form a smooth compact null hypersurface, i.e. 
a Cauchy horizon. It was shown in [23] that if the spacetime is an analytic solution of 
the vacuum field equations and if the Cauchy horizon is ruled by closed null generators 
in the sense of an S^-bundle, i.e. the null generator coincides with the generators of 
the bundle, then the spacetime necessarily has a 1-dimensional isometry group and the 
corresponding Killing field is proportional to the null generators of the Cauchy horizon 
on the horizon. The result was generalized to the smooth case in [T3] . 

Motivated by this result, Moncrief focuses his discussion of generalized Taub-NUT 
solutions on the case of spacetimes above for which the Cauchy horizon is generated in 
the sense of a Hopf bundle with U (1) isometry group. The metrics of all such spacetimes 
can be written as 



for sufficiently small t > 0, where we can assume as before that the corresponding Killing 
field is 9p^. The functions 7 and only depend on t, 9 and p2- The index a takes the 
values 1 (corresponding to the coordinate 9) and 3 (corresponding to the coordinate 
p2)- The field cjab is a symmetric 2-tensor field and the Pa a 1-form. Moreover, > is 
a constant. One assumes that all fields 7, iV, gab and (3a are smooth on (0, 5] x §^ for 
some small 5 > and have a unique smooth extension through t = 0. We express this 
as the requirement that these can be considerecl^ as smooth fields on [—5, 5] x S^. Here, 
we assume that N is strictly positive. 

In the following, we restrict to the case of Gowdy symmetry. Therefore the metric 
coefficients above are in addition independent of p2- As discussed before, the metric can 
then be assumed to be of the form Eq. ((Tj) with Eq. (ITT]) for t G (0, 5] for some 5 G (0, tt). 
Let us defin^ the function by 



^ This is not to be confused with the requirement that the metric g is smooth on [—5, S\ x S'^. Indeed, 
(7, in the form above, degenerates at t = and hence cannot be extended as a Lorentzian metric to all 
t e [-5,5]. 

^ Eqs. JTC]) and ^ imply that N{t, 0) = N{t, tt) = 1 for all t > 0. 



9 = 



e-^^-N^dt^ + gabdx^dx^) + smHe^^[k{dpi + cos^dps) + Padx""]^, 
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Then we must set 

gee — e , gep^ — 0, gp^p^ — — — — , 

(3e =0, (3^ = RoiQ~cos9), 
,7 RoN M RoN 

^ = ^' ^ = -T + ^"^- 

The next task is to identify the conditions under which these spacetimes, in 
particular the metrics g, are extendible smoothly through t = 0. We introduce new 
coordinates {t' , 9' , pi , P2') by 

t = arcsin Vf, 9 = 9', pi = p/ + ^ In t', p2 = p-l ■ (29) 

Rq 

Here, we assume that t > is sufficiently small and hence t' > is small. The quantity 
K is a constant which has, so far, not been fixed. In these new coordinates, the metric 
becomes 

/ M -Mivr2 2\ 

g = - - ' ^ dt'' + e^^d9' 



^4(1 - t')t' t 
+e-*^iV2 



2i?o^(dpi' + Qdp2')dt' + R^t' (dp/ + Qrfpa') J + ^^dp2'' 

The metric extends smoothly through t' = if the functions M, A^^ (strictly positive) 
and Q are smooth with respect to t' in a neighborhood of t' = and if (4/t^A^^ — e^^'^)/t' 
is smooth for some choice of the constant k. In this case, the field dpi^ = dp^ is a 
null generator of the surface given by t' = 0. This surface is therefore a smooth null 
hypersurface with S'^-topology and closed null generators and so is a smooth past Cauchy 
horizon. 

For the Taub solutions for example, we choose 



K = ± + m(m + \/P + ra^] 

Then the metric is extendible analytically and hence the t = 0-surface is a past Cauchy 
horizoijfl. The analogous argument applied for t = n implies the existence of a future 
Cauchy horizon. It is shown in [311 [12] that several non-equivalent extensions through 
both horizons exist. 

We will call spacetimes with all the above properties smooth Gowdy symmetric 
generalized Taub-NUT spacetimes. The name is motivated by Moncrief's before 
mentioned notion of generalized Taub-NUT spacetimes. We will restrict to Gowdy 
symmetry, but in contrast to Moncrief, we will not assume analyticity. Note that 
if an analytic spacetime as above solves Einstein's field equations in vacuum with 
vanishing cosmological constant for t > 0, then the analytically extended spacetimes 
are necessarily also solutions of the vacuum field equations. In the non-analytic smooth 
case, we do not know in general whether the extensions are vacuum solutions. We will 
not address this problem in this paper. As an interesting side-remark: Chrusciel et al. 

^ Often in this paper, we will speak sloppily of the = 0-surface" when we actually mean the t' = 0- 
surface. 
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note in [9] that there are no smooth extensions through a Cauchy horizon of S^-topology 
- solution of the field equations or not - in the polarized Gowdy case. As we mentioned 
before, however, none of the spacetimes which we consider in the following are polarized. 

3.3. Existence of smooth Gowdy symmetric generalized Taub-NUT solutions 

3.3.1. The main existence result. In this section we address the question as to whether 
smooth Gowdy symmetric generalized Taub-NUT spacetimes, which solve Einstein's 
vacuum field equations, exist. Our techniques here are based on the particular Fuchsian 
method introduced in [3 |6]; a quick summary can be found in the appendix. Similar 
existence results can be obtained by the Fuchsian techniques in [211 ESI ES]- 

In the following we call a function rotationally symmetric on if it does not 
depend on the azimuthal angle cf) in standard spherical coordinates Eq. (|6]). The Hopf 
map allows to lift any such function to a smooth f/(l) x f/(l)-invariant function on 
(or R X if the function also depends on t). 

Theorem 3.1 Let S^^: and Q* be rotationally symmetric functions in so that 



and Rq a positive constant. Then there exists a unique smooth Gowdy symmetric 
generalized Taub-NUT solution for all t G (0, vr) with 

R(t, 6) = Rq sin if: sin 9, 

satisfying the following uniform expansions at t = 0: 



Gorresponding expansions hold for all derivatives. 

Let us make a couple of comments. 

(i) This result means that we can construct a unique smooth Gowdy symmetric 
generalized Taub-NUT solution from any given asymptotic data functions 5**=,, and 
Q* subject to the conditions (130|) . We have thus obtained the same number of free 
functions as in Moncrief 's class of generalized Taub-NUT solution (after factoring 
out gauge transformations in his class). 

(ii) The Taub-NUT solutions are determined by the asymptotic data 



^^^(O) = ^^^(Tr) 



(30) 



i?(t,^)e^(*.^)=tV"(^)+0(t^), 

Q{t,e) = cose + Q,{e) sin^e + o{t^), 

M{t, 9) = S,,i9) - 2S„(0) + 21n/?o + 0{t^). 



Ro = 2lVP + m'^, 




Q. =0. 
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(iii) We shall now show that none of the solutions of Theorem 13.11 are polarized. 
Recall from Section 12.41 that a smooth Gowdy symmetric spacetime is polarized 
if Q vanishes with respect to the (Ai, A2)-parametrization of the symmetry orbits. 
Eq. f l23p yields that the spacetime is polarized if and only if 1 — (1 — Q^)e^^ = with 
respect to the (pi, p2)-paranietrization. However, for our solutions, Q is bounded 
in a neighborhood of t = and e"^^ is 0{t^) for every 6, which is therefore a 
contradiction. 

(iv) We compute the uniform limit of the quantity N (defined in Eq. (125]) ) at t = 0, 

iV(0,^)=e2(^"(^)^^"(°)). 

From that it is easy to determine that value of the constant n (defined in Eq. fl29|) ) 
which allows the solution to be extended through t = 0: 

2 

3.3.2. Equations and unknowns. In order to attempt the proof of Theorem 13. ![ let us 
make the following convenient choices. Instead of the unknown L, we shall use 

S ■.= \xi\ = L + \xiR. 

Moreover, Q is replaced by the twist potential u (defined in Section [2]). Eqs. ([1]) and 
([2]), together with f[T^ and ffT^ imply the following equations for 5* and cj, 

D'^S - t^A§25 ={l-t cot t)DS - e-'^^[{Dujf - {tuje)% (31) 

D'^uj - ADuj - t^Agaw = (1 ~tcott)Du + 2{DS - 2)Duj - 2{tdgS){tdgUj). (32) 

Note that we have added a term —ADco to the second equation for later convenience. We 
use the notation D := tdt and = tdtitdt). The operator A§2 is the Laplace operator 
of the standard metric on the unit sphere 

A§2 = dl + cot Ode + -^dl 

In our case all unknowns are independent of the azimuthal angle 0, i.e. are rotationally 
symmetric. Therefore these two equations are geometric wave equations with respect 
to the standard metric on the unit sphere as long as t G (0, tt). Indeed, these equations 
form a second-order hyperbolic Fuchsian system, see below. 

Eq. ([3]) gives us the remaining equations. On the one hand, we obtain a wave 
equation for M, 

= -Mtt + A§2M + cot t{Mt + 2St) - S'^ - e'^^w^ + 2 =: (33) 
On the other hand, we get two first-order equations 

= 4i?±,e + R{Sl + e-^^ool) - 2R±{S± + M±) =: 2i?±C±, (34) 

witli^9-|- := d0±dt. The quantities H, C+ and C_ are introduced for the later discussion. 

^ Lower indices ± of any function previously defined without an index mean the directional derivatives 
along the vector fields d± . 
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3. 3. 3. Steps of the proof of Theorem \3.1[ The main step in the proof is the construction 
of solutions of flSTl) and (1521) in a small time neighborhood of t = 0, which are compatible 
with the notion of smooth Gowdy symmetric generalized Taub-NUT solutions given in 
Section 13.21 The main local existence result is as follows. 

Proposition 3.2 Let G M, and u;**,^** G if'^^^(S^) he rotationally symmetric 
functions with k > 3. Choose constants a = (ai, 02) with 1 < ai < 2 and < 0:2 < cti- 
Then there exists a unique solution {S, u) of the Gowdy equations l[31\) and [3^] with 

S{t,d) = 2\nt + S,,{d)+wi{i.d). 
u{t, 9)=uj, + uj,,{9)t^ + W2{t, 9), 

where w := (^1,^2) G Xs^a,k for a sufficiently small 6 > 0. The functions Wi,W2 are 
rotationally symmetric for each t G (0,5]. //a;** andS^,^ are in C°° (E)"^) , thenw G Xs,a,oo- 

The spaces Xs,a,k and Xs,a,k are introduced in the appendix. The spaces are the 
standard Sobolev spaces. 

Here are some comments about this result: 

(i) The problem addressed in this proposition has many similarities to the standard 
singular initial value problem of Fuchsian equations. Hence Theorem 3.6 of [B] 
should apply directly. However, in contrast to the standard singular initial value 
problem, it is crucial here to assume that = constant. Only if this is the case, 
are conditions (1) and (2) of that theorem satisfied. Still, condition (3) fails. We 
find that condition (3) has been derived in 0| from Proposition 3.5 for non-constant 
asymptotic data. Given that w^, = constant, however, condition (3) weakens and it 
is straightforward to check that these assumptions are sufficient. 

(ii) We stress that the singular initial value problem considered here is significantly 
different than a standard Cauchy problem. Since the equations are singular at 
t = 0, but of Fuchsian type, we obtain the fourth power of t in the leading-order 
term of u and the logarithm in the expansion of S on the one hand, and we do not 
obtain the full expected number of free functions on the other hand. However, as 
mentioned before, the number of free functions here is the same as in Moncrief 's 
results. 

(iii) We must assume > 3 (as opposed to k > 2 in [6]). Because of the Sobolev 
inequalities, the condition A; > 3 is required when we check condition (2) of 
Theorem 3.6 of [6] for the 2-dimensional spatial manifold here. It is crucial 
that Eqs. fl?I]) and fl52]) are geometric wave equations in order to derive the energy 
estimates in the same way as in [7',^. 

(iv) The last statement in the theorem means that if w** and S"** are smooth rotationally 
symmetric asymptotic data functions on then the corresponding solution is 
smooth on the time interval (0,5]. Furthermore, w and all of its derivatives decay 
with the same rate in t at each spatial point. 
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(v) The constants ai and a2 control the decay of the remainders wi and W2 at t = 0. 
Roughly speaking, the statement {wi,W2) G Xs^a,k means that Wi = 0(t"i) and 
W2 = 0(t''^"^) at t = 0. We will argue below that in fact ai = 02 = 2 in the 
smooth case. The fact that the statement about these constants in Proposition l3.2l is 
weaker than this is related to the following property of our general Fuchsian theory. 
Roughly speaking, the theory assumes that the solutions might have additional 
logarithmic terms in t. Those terms would show up in particular if the exponents 
of t in the leading-order terms at t = are spatially dependent. This is generically 
the case for many Fuchsian problems, however, not for us here. 

The proof of Proposition 13.21 follows roughly the same steps as in the Gowdy case with 
spatial T^-topology discussed in [HI |6]. 

Let us restrict to the smooth case u^^, S^^ G C°°(S^) from now on. We see easily 
from Eqs. f l3T|) . f l32|) and the expansions in Proposition 13. 2l that e*^ and u can be extended 
as smooth solutions through t = such that both are even in t. Hence, we can assume 
in all of what follows that both are solutions of the equations on, say, [—5, 6] x S^. 

We can compute Q by means of Eqs. 

deQ = -Re-^^dtu, dtQ = -Re'^^deco. (35) 

The integrability condition for these equations follows from the fact that w is a solution 
of Eq. 032]) • Let us suppose that a solution of Proposition 13.21 is given. Then all terms 
in Eqs. (135|) are smooth on [—6, 6] x Hence these equations determine Q uniquely 
up to a constant and Q is a smooth function on [—6, 6] x S^. From 

deQ{0,9) = -ARoUj,,{e)e-''^**^'^ sin^, 

it follows that 

Q{0,9) = Qo-4Ro / a;„(x)e-2^"(")sinxda;, 
Jo 

for some constant Qq. In order to guarantee the smoothness condition Eq. f|T7|l . namely 
Q(t, 0) = 1, Q{t, tt) = —1 for all t > 0, we must choose Qq = 1 and the asymptotic data 
must satisfy 

w**(a;)e~^'^**''^'' sinx dx = . (36) 

2i?o 

Then it follows that (5(0,0) = 1, (5(0, tt) = —1. Hence from the second of Eqs. fl35l) we 
conclude that Q{t,0) = 1, Q{t,7c) = —1 for all t G [—5,5]. This is summarized in the 
following lemma. 

Lemma 3.3 Let Rq > be a constant and w** and 5*^,* be rotationally symmetric 
functions in (7°°(§^) satisfying Eq. (EQj- Suppose that {S,ijj) is the corresponding 
solution of / Tfi]) and [3^) according to Proposition \3.^ Then there is a unique smooth 
solution Q of Eqs. i\35\} on [—6, 6] x §^ which is rotationally invariant at each t and 
Q(t,0) = 1, Q(t,7r) = -1. 
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It follows easily that Q has the expansion 

Q{t,e) = l-ARo / a;„(x)e-2^"(")sina;dx + -i?ow:*We-''^"('^sin^t2^..., 
Jo ^ 

in t at t = and that Q is an even function in t. We can introduce a smooth rotationally 

symmetric function Q^{6) such that 

Q{t,e) = cosO + Q^e) sin^ e + 0{t'^). 

Hence, instead of prescribing the function as one of the asymptotic data functions 
in Proposition 13. 2[ we can equivalently prescribe the function and then set 

1 - d,^Q, (6) sin 6 - 2Q,{e) cos 6 



4i?n 



In the next step of the proof of Theorem 13 .11 we construct the function M from the 
constraint equations flM|) for given solutions (5, u) of Proposition 13.21 The result is as 
follows. 

Proposition 3.4 Let Rq > be a constant and w** and S*** be rotationally symmetric 
functions in C°°(S^) satisfying 

Suppose that {S, u) is the corresponding solution of i\3l\} and [3^) according to 
Proposition \3.2[ Then there is a unique smooth solution M of Eqs. [37\ ) on [—6,6] x 
which is rotationally invariant at each t and satisfies the smoothness condition e*^*^*'^) = 
(/^g sin^ t)e-^(*'^) for all {t,9) e [-6,6] x {0,7r}. Moreover, M is an even function oft. 
The expansion of M at t = in t is 

M{t, 6) = S,,{e) - 25„(0) + 2 lni?o + 0{t^). 
The proof of this is maybe the trickiest part of this section. Let us define the functions 

/i± ■.= ^R±,e + R{Sl + e-^^ojl). 
Then the constraint equations can be written as 

Since R± = Rosm(t ± 9), the terms on the right-hand sides are potentially dangerous 
along the diagonals of the Gowdy square where R± = 0. Moreover, the functions fi± 
are not bounded everywhere in the limit t — ?► 0. We proceed as follows. Instead of 
interpreting the functions in Eqs. f l37|) as rotationally symmetric functions on M x 
we consider them for the moment as functions on M x M where we extend 6 to M, 27r- 
periodically. Smooth rotationally symmetric functions on M x will then be identified 
with functions which are even periodic in 6. Since S is such a smooth function, it follows 
that St is even and Sg is odd in this sense. Since we want to construct a smooth solution 
M of Eqs. §7i), it follows that 



AR, AR_ r 4R. AR_ 
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are even and odd functions of 6, respectively. Indeed this is true and it follows directly 
from their definition. The term 2/t is introduced for later convenience. This renders 
Eqs. (EZD into 

Mt = ~{St - 2/t) + F"^"^, Me = -Se + (39) 

The fact that {S, u) solves Eqs. (131 p and fl32p implies the following linear symmetric 
hyperbolic system 

with the even source function 

Note that this function extends smoothly through t = and hence the system Eq. (1401) 
has a well-posed Cauchy problem with data on t = 0. We call those data functions 
F^^'^^{9) and F"'^'^{9), and hence must suppose that these are even and odd, respectively, 
periodic smooth functions. Let us write FJ™"^ and F^'^'^ for the smooth solution of 
Eq. (HOj) for any given such data and in order to distinguish them from the 

functions F^^^"^ and F°'^'^ determined by Eq. (155]) from the given functions S and u. It 
is straightforward to solve Eqs. (1391) for M when _p*^^^° is substituted by FJ^'^" and F°'^'^ 
by F^'^'^. Then the integrability condition is satisfied as a consequence of Eq. fHU]) . and 
it follows 

M{t,9) = M,-{Sit,9)-2lnt)+ [ F°'^'^{x)dx+ [ F;™'^(r, 0) dr, (41) 

Jo Jo 

for some constant M*. Therefore, M can be considered as a smooth rotationally 
symmetric function on [—6, 6] x S^. It is an even function in t if and only if F^"^™ = 0. 

Now we must study the conditions for which .F'^™^ = ^even jpodd ^ ^odd -^^ 
order to give a meaning to Eq. (14T|) . We find that under the conditions for S and u 
of Proposition 13. 2[ the function .F^™" converges pointwise to at t = 0, while 
goes to 25*^^, for every 9 G [0,7r]. Indeed, F"™'^ and F"'^'^ can be extended uniquely 
as continuous functions to all the points given by t = and 9 G (0,7r). Let us hence 
choose F^^^"^ = and F^'^'^ = 2S'^^. On the interior of the domain of dependence for 
Eqs. (HOD of the subinterval (0, vr) of the t = 0-surface, it follows that F^™'' = Fp"" 
and F'^'^'^ = Fj'^'^ because both sets of functions are smooth and satisfy Eq. fHOl) with 
the same data at t = 0. Let us call this domain Q. Recall, however, that F*^^^" and 
^odd g^j^g (jggned in those points of dfl where R±. = 0. However, we can extend 
the functions F'^'^™ and F"'^'^ to those points using the continuous values of FJ™"^ and 
Fj'^'^ there. From the definition of F'^™'^ and F°'^'^ involving only quotients of smooth 
functions, which become zero simultaneously at the same points, we obtain that these 
extensions of F''^'^" and F°'^'^ are continuous on the whole domain, in particular through 
dQ, possibly except for the points (0,0) and (0,7r). The continuity in the points (0,0) 
and (0, vr) follows directly from the asymptotic behavior of 5* and u there. Therefore, 
^even jpodd uniformly continuous functions on the whole domain. Now, both 
sets of functions F^™'^, F°'^'^, and Fp°, Ff'^, are solutions of Eqs. (gO]) with the same 
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smooth data at t = 0. It is a standard fact about symmetric hyperbolic systems that 
corresponding solutions are smooth and that there can only exist one smooth solution 
for given smooth data. Therefore, it follows F'^'"^^ = i^j™"^ and F°'^'^ = F^'^'^ everywhere. 
So, Eq. ( HTj) represents the actual solution of Eq. (15^ . which takes the form 

M{t,e) = M, + 2S,,{e)-2S,.{0)-{S{t,e)-2lnt)+ [ F'=™'^(r,e)dr.(42) 

Jo 

Now let us consider the smoothness condition e^*^*'^-* = (R^ sin^ t)e~^^^'^^ for all 
t G (0, 6] and 6* = 0, vr. We compute 

dt{M-2\nRo + S-2\n sin t) = F^™°|e=o,. + ^ - 2 cot t. 

Evaluating .F'^™'^ from the definition Eq. f l38p at 6' = 0,7r (where in particular R = 0), 
we conclude that 

dt{M -2\nRo + S -2lnsmt)\g^Q^^ = 0. 

Therefore, the smoothness condition for M is satisfied for all times if and only if it is 
satisfied at t = 0. From Eq. (H2l) . we see that this is the case provided M* = 21n_Ro and 
5'=i,^=(0) = S^^{7i). This concludes the proof of the last proposition. 

Now we define Ci := (C+ + C_)/2 and C2 := {€+ - C_)/2 from Eqs. ([21, and 
is defined in Eqs. fl33|) . The system fl3T|) and fl32|) implies the subsidiary system 

dtCi - d^C2 = 0, dtC2 - d^Ci = H + cottC2 + cot 9 Ci. 

Since M is a solution of Eqs. (l37j) . it follows that the quantities Ci and C2 vanish 
identically. Then H must also be zero. We conclude that our solutions indeed solve the 
full set of Einstein's field equations. 

All the so far constructed functions can now be lifted to smooth functions on 
which are invariant along dp-^ and dp^ at each t = 0. All these functions are smooth 
functions on [—5,5] x S'^ with respect to the coordinates (t' ,9' , p[, p'2) for some small 
S > 0. We have thus obtained smooth Gowdy symmetric generalized Taub-NUT 
solutions in a (possibly small) time interval < t < 6. The global existence theorem 
of Chrusciel [1^ can now be used to extend the spacetimes to the whole time t interval 
(0, tt) as smooth globally hyperbolic Gowdy solutions. 

4. The linear problem and global-in-time properties 

We have seen above that for given smooth asymptotic data at t = (e.g. the values 
of S*** and u;**) a smooth Gowdy-symmetric generalized Taub-NUT solution exists in a 
vicinity of t = 0. Moreover, using Chrusciel's theorem (theorem 6.3 in [ID]), we see that 
this solution can even be extended smoothly to the whole time interval (0, vr). However, 
the surface t = ir itself is expected to contain either singularities or Cauchy horizons. 
On the other hand, Chrusciel's result also allows the case that the t = vr-surface is 
regular, but just the coordinates break down there. It is the purpose of the following 
considerations to find out what happens at t = tt. In particular, we construct explicitly 
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the metric potentials at this boundary (as well as on the axes 6 = 0, tt) in terms of 
the asymptotic data. For that purpose, we apply the so-called soliton methods, which 
were used in [19] for the investigation of x Gowdy spacetimeJ^. to the present 
§^-symmetric case. 

In all of what follows we make the same hypotheses as in Theorem 13. 1[ These 
assumptions are consistent with those listed in Section 13. 2[ and hence we consider 
"smooth Gowdy symmetric generalized Taub-NUT solutions" . 



4.I. Einstein's field equations and the Ernst formulation 

The first important step for the following considerations is the introduction of the 
complex Ernst formulation of the Einstein equations which will be described in this 
subsection. 

Again we start from the metric 

g = e^^i-dt"^ + d^^) + i?osintsin^ [e^(dpi + Qdpa)^ + e'^dpl] (43) 

in the Killing basis {dp-^, Spg}. Here, we express L in terms of a metric potential u via 

e^ = ^e«. (44) 
sm (J 

In this way, we arrive at 

g = e*^(-dt^ + d^^) + Ro [sin2te"(dpi + Qdpa)^ + sin2^e~"dp^] . (45) 
Note that u is related to the quantity S (defined in Sec. I3.3.2p via 

u{t, 9) = S{t, 6) - ln(/2o) - 2 In sin t, 

i.e. the singularity of 5* at t = (S* behaves as 2 In t for t — 0, see Prop. 13. 2p is removed 
by subtracting the term 2 In sin t. 

Now we reformulate the Einstein equations as equations for m, Q and M. We obtain 
two second-order equations for the metric potentials u and Q, 

- d^u - cot t dtu + dlu + cot d deu + e^" — ^ \{dtQf - [deQf] +2 = 0, (46) 

sin 

- - 3 cot t dtQ + dlQ - cot d deQ - 2[{dtu){dtQ) - {deu){deQ)] = (47) 

and two first-order equations for M, 

1 sin'^t r 

(cos^t - cos^9)dtM = -e^"^^ costsm9[{dtQf + {deQf] - 2 sin t cos 9 (dtQ) (deQ) 



2 sin^ 

H — sin t sin 9 
2 

— (2 cos^t cos^6' — cos^t — cos^6') dtU 



cos t sin 9[{dtu) + (deu) ] — 2 sin t cos 9 {dtu){deu) 



— 2 sin t cost sin 9 cos 9 {deu + tan^^), (48) 

The methods described |19) have also been apphed for studying the interior region of axisymmetric 
and stationary black holes with surrounding matter, see P1I51IT5]. 
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1 sin'^f r 

(cosH - cos^e)deM = - -e^"— — sint cos e[idtQf + (deQf] - 2 cos t sin 9 {dtQ){deQ) 



2 sin^ 

sin t sin 6 

2 



dta = — sin if: sin 9 deb, dga = — sin t sin 9 dtb (52) 



sintcos^^[(9tu)^ + (deu)^] — 2 cos t sin 9 {dtu){deu 

— 2 sin t cos t sin 9 cos 9{dtu — tan t) 

- (2 cos^t cos^^ - cos^t - cos^^) deu. (49) 

Since M does not appear in f l46p and 0471) and since we assume the genericity condition 
of Chrusciel, these equations may be solved in a first step. Afterwards, f HHjl and fH9|l 
can be used to calculate M via a line integral. Note that the integrability condition 
dtdgM = dgdfM of the system fHS]) . fH^ is satisfied as a consequence of fH^ . fH7|) . 
Hence, M does not depend on the path of integration. 

It turns out that the two Einstein equations (H^ . (jUj) are equivalent to a single 
complex equation, namely to the Ernst equation 

/ {-d'^S - cot t dtS + dl£ + cot 9 deS) = -{dtS^ + {de£f (50) 

for the complex Ernst potential S = f + ib. Here, the real part / of £^ is defined in terms 
of the Killing vector dp^ by 

/ := -^9idp,,dp,) = Q^e" sin^t + e"" sin^^ (51) 

and the imaginary part b is given by 

— sin t sin 9 dgb, dga = — 
r J 

with 

„:=9§44 = «e«sm^(. (53) 

Note that for smooth functions u and Q the Ernst potential S is also smooth: For the 
real part /, smoothness is clear from definition (IHT]) . In the case of the imaginary part b, 
it can be shown by solving the two equations in fl52l) for dtb and dgb and replacing a and 
/ via flS^ and fl3T]) . The resulting expressions for dfb and dgb in terms of Q and u (and 
their first order derivatives) turn out to be smooth functions, if we use the fact that Q 
behaves as given in 0171) . Hence, integration will lead to a smooth function b. Therefore, 
we can conclude from the previous local existence results and Chrusciel's global existence 
theorem that for any given set of asymptotic data (as described in Theorem 13. ip the 
corresponding Ernst potential £^ is a smooth complex function on (— T, it) x §^ for some 
T > 0. In the following we investigate under which conditions S can be extended 
smoothly to the boundary t = n and beyond. Note that our assumptions imply that 
/ > holds in the entire Gowdy square with the exception of the points A and B and 
with the possible exception of the future boundary "Hf (see Fig. [1]) which is important 
since we will divide by / in some of the following formulae. 

Once we have obtained an Ernst potential £^ as a solution to the Ernst equation 
f l50p . we can calculate the corresponding metric potentials from it. It turns out that 
the integrability condition dfdga = dgdta of fl52|) is satisfied as a consequence of the 
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Figure 1. Wc integrate the LP along the boundaries of the Gowdy square (dashed 
path) in order to investigate for which asymptotic data the solution can be regularly 
extended up to the future boundary "Hf {t — tt). 



Ernst equation. Therefore, a may be calculated via line integration from S. The metric 
potentials u and Q can then be obtained from a and /. With (15T|) and (153|) we find 

e"-^ + ^. Q- „./:\... (54) 

sm t J j'^a'^ + sm ism U 

Finally, M may be calculated using ( HSj) and (HH]), as mentioned earlier. 
As an example, we give the Ernst potential for the Taub solution: 
11 X 

f=— sinH cos^e + — sin^e, (55) 
jC 21 

1 



cost(cos^t — 3)Vm'^ + P — 2m cos^6' + cost (56) 



with X := (1 + cosH)Vm'^ + P + 2m cost. (Here we have set an arbitrary additive 
integration constant in b to zero.) 

4-2. The linear problem 

Interestingly, the Ernst equation (150|) belongs to a remarkable class of nonlinear partial 
differential equations for which an associated linear problem (LP) exists which is 
equivalent to the nonlinear equation via its integrability condition. For applications 
of this LP in the context of axisymmetric and stationary spacetimes we refer the reader 
to, e.g., ED]- In the Gowdy setting, we use the LP in the form ES], which reads 
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\ / 

A, r \ A, 



By \ 1/ By 

Ay X\ Ay 



(57) 



where the pseudopotential $ = y, K) is a 2 x 2 matrix depending on the coordinates 

x = cos(t-6'), y = cos{t + e) (58) 

as well as on the spectral parameter K E C The remaining ingredients of the LP are 
the function A, 



X{x,y,K) :-- 



K 



y 



(59) 



A, 



x,y. 



(60) 



K-x' 

and the matrix elements A^, Ay, B^ and By, defined in terms of the Ernst potential as 

djS _ djS 

2/ ' 2/ ' 

Due to the two possible signs of the square root in (159|) . A : C — t- C, K h-j- A 
describes, for fixed values x, y, a mapping from a two-sheeted Riemann surface (K- 
plane) onto the complex A-plane. The two if-sheets are connected at the branch points 

Ki=x (A = cx)), K2 = y (A = 0). (61) 

In general, the pseudopotential $ will take on different values on the two fT-sheets. 
Only at the branch points it has to be unique, since both Riemannian sheets coincide 
there. We will see below that this observation plays an important role for the calculation 
of the Ernst potential from the solution of the LP. 

As already mentioned, the integrability condition dxdy^ = dydx^ of fl57|) is 
equivalent to the Ernst equation flSU]) . Hence, the Ernst equation is a consequence 
of the LP and, on the other hand, for a given potential £^ as a solution to the Ernst 
equation, the matrix $ does not depend on the path of integration. 

Finally, we note that for any solution # to the LP (lFr|) . the product ^C{K), where 
C{K) is an arbitrary 2x2 matrix, is also a solution (corresponding to the same Ernst 
potential). As shown by Neugebauer [2S], it is always possible to choose C{K) in such 

The formal relation between our coordinates (describing Gowdy spacetimes with two spacelike Killing 
vectors) and the Weyl-Lewis-Papapetrou coordinates {pX,(p,i) as used by Neugebauer (describing 
axisymmetric and stationary spacetimes with one spacelike and one timelike Killing field) is given by 
p = ii?o sin < sin 6*, C = i?o cost cos 0, ip = pi, t = p2. 
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a way that the transformed pseudopotential takes the form 



i^i{x,y,K) i)f{x,y,K) 



^<{x,y,K) 



i^iix,y,K) ^ljf{x,y,K) 



(62) 



1 
-1 



^>{x,y,K) 




where the superscripts ">" or "<" indicate whether the functions are evaluated on the 
"upper" (A = 1 for K = oo) or "lower" (A = —1 for K = oo) i^-sheet. Hence, $ is 
completely determined by the values of two functions ipi and ip2 on both i^-sheets. In 
all of what follows we assume that we have already achieved this form of 



4-3. Solution of the linear problem 

4-3.1. Coordinate transformation. In the following we intend to integrate the LP along 
the boundaries of the Gowdy square. For that purpose, it turns out to be useful to study 
the situation not only in the coordinate system E, corresponding to the Killing basis 
{c?pj,9p2}, but also in a coordinate frame S, 

S : i=t, 9 = 9, pi = pi + qp2, p2 = p2 (63) 

with q = constant. According to f|T8|) - fl20|) and (jS]), the transformed metric potentials 
are 



Ro — Ro, 



u = u, 



Q = Q-q, 



(64) 



i.e. only Q is changed by subtracting a constant. In particular, we will choose the two 
systems S with g = lorg = — 1, in which QUi = or QUa = holds, respectively. 

Since the coordinate transformation fl63l) is merely a change of the Killing basis, 
the Ernst equation fl50l) retains its form in S. This implies the existence of a LP fl57j) 
for a pseudopotential $ in this frame. As shown by Neugebauer [271 EH], the matrices 
$ and # are connected by the transformation 



with 








x) 



1 

-A 



A 



a ± y sin t sin 9 



(65) 



(66) 



where all quantities on the right hand side of Eq. fl65|) belong to the original frame S. 
As we will see, this transformation becomes particularly simple at the boundaries of the 
Gowdy square for our choices g = ±1. 



Smooth Gowdy symmetric generalized Taub-NUT solutions 



25 



4-3.2. The LP on Hp, Ai and A2. From our previous discussions, namely from the 
local investigation of the singular initial value problem for the Einstein equations with 
Fuchsian methods (see Sec. 13.31) and from Chrusciel's global existence theorem, we 
know that for any smooth set of asymptotic data on "Hp a corresponding smooth Gowdy 
symmetric generalized Taub-NUT solutions exists. Moreover, this solution is smooth 
both on the axes of symmetry Ai, A2 and in the interior of the Gowdy square. The goal 
of this subsection is to find explicit expressions for the values of the solution on Ai and 
A2, which are determined by the data on "Hp. Afterwards we may study the behavior 
as t — )■ TT and investigate whether a continuation of the solution to "Hf is possible. 

Along the entire integration path, we have x = y and therefore A = ±1, cf. fl59|) . 
However, it suffices to study the case A = 1 alone, since the solution on the Riemannian 
sheet with A = — 1 can easily be obtained from the solution with A = 1 using (l62l) . 

For X = y and A = 1, the LP (1571) reduces to the ODE 




dx£ dx£ 
dx£ dx£ 



$ (67) 



with the general solutior 

$ = EC(if), E := I ; ^ I (6^ 




in terms of the Ernst potential on the boundary, where the 2x2 matrix C is a in- 
dependent "integration constant". The solutions on all parts of the integration path 
have the form fIBS]) . but with different integration constants: 



« = : # = EC, C = X' . (69) 




« = 0: # = ED. D=| ^^j. (70) 

« = - * = ED, DH (71) 

A further simplification can be achieved by normalizing $ at t = via 

t = 0: ^<=^< = i;(K), (72) 

where ip{K) is an arbitrary gauge function (which will later be specified in such a way 

that the LP has a regular solution, see Eq. (183|) below). This is possible since the form 
(1^2]) of $ is invariant under the transformation 



' P{K) a{K) ' ^ ' 

By plugging the solution (1551) into ([57)) . we see that the matrix on the right hand side of (j67p is 
proportional to /, i.e. the factor 1// is cancelled and the solution extends smoothly over points with 
/ = 0. However, it follows from (l5T|) that / does not vanish on "Hp, Ai, A2 with the exception of the 
corners B of the Gowdy square (as well as C and Z?, provided Q^e" is bounded for t — tt). 
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The two degrees of freedom a, (3 can be used to achieve the two conditions in ( 1721) . As 
a consequence, we obtain 

C, = 0, C^ = 4j (74) 

in this gauge. 

From ( l69|) - (m]) we can now calculate the solution of the LP in the frame S (cf. 
f l63p ) using the transformation formula fl65l) . It follows from (1531) that a takes on the 
boundary valued 

"Hp : a = 0, Ai : a = — = 1, A2 ■ ^ = = i^^) 

Plugging this into fl65|) . we obtain for A = 1 



9 = 0: ^ = +2i{K -x)\ in S with g = 1, (77) 

9 = 71 :^ = -2i{K - x) I ^~ ^~ I in S with q = -1. (78) 

-Di —D^ J 

As we will see below, Ci{K) and C2{K) are determined completely by the data at 
t = 0. Now we intend to express the components of the matrices D and D in terms of Ci, 
€2- For that purpose, we use that $ has to be continuous at the corners A and B of the 
Gowdy square (see Fig. [1]). This condition leads to an algebraic system of 4 equations 
which, however, is not sufficient to calculate the 8 unknowns Di, . . . , D4, Di, . . . , -D4. 
This is the reason for introducing the coordinate frame S. From the requirement that 
also # (in E with g = 1) is continuous at A and $ (in E with g = — 1) is continuous at 
B we find another 4 algebraic equation^. In this way we obtain an algebraic system of 
8 equations for the 8 unknowns with the following solution for the matrices 'D{K) and 
t>{K) in terms of C{K): 

n ^ ^ACi + iCs ^ „ ihA{hACi + iC2) 

— '^l i^TTTT ^2 — ^2 



2{K-1) ' 2{K-1) 



2(K-1)' " ^ V 2(K-1) 

-t^l — H ^777 7T— , J^2 — O2 + 



2(ir + l) ' 2{K+1) 

Note that a is automatically discontinuous at the points A, B, C, D as a consequence of the definition 
(|53l) . In contrast, a is a smooth function in the remaining part of the Gowdy square. 

The reason why just the usage of a different coordinate system can lead to independent algebraic 
equations is the following. The boundary values of the quantity a enter the transformation law (j65p . 
Therefore, the additional algebraic equations found in E ensure that a indeed takes on the boundary 
values (I75p (and, as a consequence, also Q takes on the correct boundary values Q = 1 on Ai and 
Q = — 1 on A2)- From (j52p alone it would only follow that a is constant on the boundaries without 
specification of these constants. 
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In the following subsection we utilize these results in order to determine the Ernst 
potential and the metric potentials on Ai and A2 in terms of the data on 7{p. 

4-3.3. Ernst potential on Ai and A2. From the solution of the LP obtained in the 
previous subsection we may, as a first step, calculate the Ernst potential on Ai and A2 
in terms of the initial potential on "Hp. To this end, we start by expressing Ci and C2 
in terms of S on Tip. 

As mentioned in Sec. \A.2\ the mapping K X in f l59|l defines a two-sheeted 
Riemannian i^'-surface. At the branch points Ki and K2, where both sheets are 
connected, any function of K has to be unique, i.e. the values on the upper and lower 
sheet have to be the same. For t-^-values on the boundaries of the Gowdy square, we 
have confluent branch points, i.e. Ki = K2 = x. The uniqueness oi ^ aX K = Ki = K2 
leads to the conditions (see fl62l)') 



7/p,A,^2: ^f=Vf and ip2 = '^2 for K = x. (79) 
In particular, on 7{p we obtain the two equations 

^pCi + C2 = ^, ^pCi - C2 = ^ (80) 
with the solution 

C,(a;) = = C2{x) = ^eMuSM^ = ^^p(^) ^.(^)^ (81) 

Sp{x)+Sp{x) fp{x)' £p{x)+£p{x) fp{x) 

where 

£p{x) = fp{x) + i&p(x) = £(t = 0,6 = arccosx). (82) 

Now we can suggest a possible choice for the gauge function ^p, which was introduced 
in dZg). If we set 

ij{K) = {K' - If, (83) 

then the solution $ (as well as $) is regular, because compensates for the poles that 
the matrices D and D would otherwise have aX K = ±1. Note that Ci and C2 are also 
regular because fp{x) = e~"(l — x'^), cf. flSTj) . But of course, as we will see below, the 
Ernst potential on Ai, A2 and Tif is independent of this gauge choice. 

With these expressions for Ci and C2 together with the solution of the LP from 
the previous section, we can also evaluate the condition (!79|) on Ai and A2 to obtain 
explicit formulae for the Ernst potential there. The result that we find independently 
of the particular choice for the gauge function ip is 

A: S.i.) -.^ea^ a^ccos... ^ 0) ^ -,^f;ltl)y 

^ / X ^/ / X ^ N \\bB — 2(1 + x)]£p{x) + b\ , , 

A2 : £2{x) := £{t = arccos -X , 6 = 7,)= \ ^ . 85 

£p[x) - i[bB + 2(1 + x)\ 
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From the latter equations we can conclude that Si and 82 are smooth functions of 
X. To see this, recall that we assume smooth data at t = 0. For smooth initial functions 
u{0,9) and Q{0,9), the initial Ernst potential £p will also be smooth, cf. Eqs. (186|) . (187|) 
below. As a consequence, the numerators and denominators of the fractions in f l84p . 
f l85|l are also smooth and an irregularity in the Ernst potentials could only occur if the 
denominators became zero for some x G [—1,1]. However, it follows from f l86|) . fl87|) 
below together with Eq. f lT7|) that the only zeros are at x = 1 or at x = — 1. Moreover, 
these equations show that the numerators have zeros of at least the same multiplicity at 
these x-values. Hence, the zeros in the numerators and denominators cancel each other 
out and the fractions are smooth functions of x for all a; G [—1,1]. The only exceptional 
cases occur for asymptotic data with 6^ = 6^ + 4 or 6^ = — 4. In the first case, £1 
diverges at point C and in the second case £2 diverges at D, cf. Fig. [TJ 



4.3.4- Metric potentials on Ai and A2. In the previous subsection we have provided 
explicit formulae for the Ernst potential on the axes Ai and A2 in terms of the initial 
potential on Tip. Now we will see how the metric potentials u, Q and M can be obtained 
from the Ernst potential on these boundaries. 

We assume that asymptotic data u{0,9) and Q{0,9) (or, equivalently, S'^,*(6') and 
and a constant Rq > are given. From these data, we may calculate the initial 
Ernst potential £p = fp + ibp. The real part can be obtained from (IHTj) . 

/p(^) = e-"(°'^) sin'e, (86) 

and the imaginary part can be calculated by integrating the first equation in fl52|) with 
respect to 9, using f l5^ . We obtain 

bp{9) = bA + 2 [ Q{0,9')sm9'd9', (87) 
Jo 

where 6a = b{0, 0) is an arbitrary integration constant. 

From £p we may calculate £1 and £2 via flM|) . (15^ . Afterwards, we can use these 
results to determine the potentials u, Q and M on Ai and A2. Using again fl5T]) - fl53|l 
together with ([12]), ([IT]) and (gl]), we find 



Ai: e"^*'"^ = ^, ,M,t,o) = i::^^ g(t,0) = l, 

A2: e"(*-) = ^, ,Mm = ^^^ g(t,vr) = -l. 
sm^t /2(t) 



(89) 



4.4- Situation on 



So far we have seen that we can prescribe arbitrary smooth data at t = and will 
always find smooth potentials for t < tt as solution to the field equations. In particular, 
we have derived explicit formulae for the Ernst potential and the metric potentials on 
the axes Ai and A2. It remains to study under which conditions the solution can even 
be extended smoothly to the future boundary Tif. In order to answer this question, we 
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tentatively solve the LP on T-Lf and investigate whether this solution can be attached 
continuously to the solutions on Ai and A2- 

The LP on "Hf reduces to the same ODE as on the other boundaries of the Gowdy 
square, namely to Eq. (1^ . We write the solution in S as 



t = n: $ = EC, C = { ] . (90) 




In order to obtain the solution in the coordinate frame S too, we need to calculate the 
quantity a on "Hf so we can apply the transformation formula (16 5p . It follows from 
that, if the metric potentials u and Q remain bounded for t — )■ vr, a = then holds on 
Tif (provided / does not vanish on Tif with exception of the boundary points C, D). 
However, it is not yet clear how u and Q behave for t — )■ vr. Therefore, so far we can 
only say that a is constant on "Hf, cf. (152|) . 

t = TT : a = aQ = constant. (91) 

Using (1651) . we find therefore in S 

t = n: inSwithg = ±l, (92) 

Now we can investigate whether $ in fl90p and $ in (1921) can be attached 
continuously to the corresponding solutions on Ai and A2- This question is equivalent 
to the solvability of an algebraic system of 8 equations. It turns out that this system 
can be solved if and only if the initial parameters Ba and bs satisfy 

bB^bA + 4, and be ^ bA - 4. (93) 

(It was already discussed at the end of Sec. I4.3.3l that the Ernst potential diverges at C 
or D if one of these conditions is violated.) The algebraic equations then fix the values 
of be, bo and ao in terms of the initial quantities &^ and 6^, 



, ^ -^bA + bBibA-be) 
bA-bB-A 



(95) 



8{bB - bA) 

«o = 777—77 j-T^- 96 

16 + {bs - bAY 

With these results we can calculate the Ernst potential on "Hf. With the same 
considerations as in Sec. 14.3.31 we obtain 

:= £{t = .,e = arccos(-x)) = ^ ^^(,) , (97) 



where 



ai = - i [{{bA - bBf + 16) - 2(6^ - &i?)(6A + - 4)a; 

+ {bA - bB? + 8(6a + bB- 2)] (98) 
a2 = 4(6^ - bB){bAbB - 2b a - 26b)x - %{b\ + 6|) (99) 
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61 = 4 [(6a - bB)x - 4] (100) 

62 = - i [{{bA - bsf + 16) + 2{bA - bB){bA + bB- 

+ (bA-bBf -8ibA + bB + 2)]. (101) 

Similarly to the discussion in Sec. 14.3.31 we find that S{ is a smooth function on the 
entire boundary "Hf (with our assumption (193|) ). 

As already mentioned, the auxiliary quantity a would satisfy the boundary 
condition a = on "Hf if the metric potentials u and Q were bounded for t ^ it. 
From fl96l) we can read off that this is only the case if the initial parameters 6a and 6^ 
satisfy the condition 

6a = 6b, (102) 

which can also be expressed in terms of the metric potential Q, 

Q{0,e)sm9de = 0, (103) 





cf. fISTI) . In the following subsections, we study separately the cases 6a = 6^ and 
6a 7^ bB. 

4.4- i- Initial data with bA = bB- Such data lead to a solution of the field equations 
with a = on Hf. As a consequence, the metric potentials u and Q are regular at Tif. 
The formula (197|) for S{ simplifies in this case to 
nj . , i{bA-l + x^)Sp{x) + b\ 

8,{x)--AbA + l-x-) 
and in terms of this Ernst potential, the metric potentials are given by 

Hi. e - ^^^^y e -Roe ^^^^^ Q[n,U) - (105j 

where /f = ^Sf, 6f = '^Sf. Here we have used that M — m is constant on T-Lf as a 
discussion of Eq. fj^9|) in the limit t — )■ tt reveals. 

It follows from these results that is a regular Cauchy horizon, generated by the 
Killing vector dp^ (just like the past horizon "Hp). To see this, we can use a modification 
of the transformation f l29p to regular coordinates in a vicinity of this boundary, 

7T — t = arcsin y/f, 6 = 6', pi = pi + — Int', p2 = P2 ■ (106) 

As a consequence of fllOSp . the constant k can always be chosen such that the metric 
is regular in terms of t', 0' , p\, p^- Moreover, [dp^.dp^ = R^e^ sinH tends to zero for 
t ^ 71, i.e. "Hf is indeed a regular null hypersurface and therefore a Cauchy horizon. 

4.4-S. Initial data with bA 7^ 65- Now we study the case bA ^ bB and assume that in 
addition 6^ 7^ 6a ± 4 holds. In this case, the auxiliary quantity a tends to 7^ as 
given in (196|) for t — vr. As a consequence of (153|) . we see that at least one of the metric 
potentials Q and u cannot be bounded in this limit. Indeed, we can read off from 0541) 
that e" diverges as 1/ sin^t for t — >• vr for all 9 G (0, vr). 
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However, it turns out that this divergence is only a pecuharity of our special choice 
of metric potentials. A better quantity for discussing regularity is the Ernst potential 
S which is defined invariantly in terms of the Killing vectors. And indeed, the Ernst 
potential also remains regular in the entire Gowdy square for 7^ bs- Moreover, from 
£ one can calculate the Kretschmann scalar RijkiR^^^'' on Tii and find that it remains 
bounded — with the exception of the earher discussed special cases 6^ = 6^ ± 4 which 
we have excluded here. (For bs — bA + 4, the Kretschmann scalar on Ai behaves as 
1/(7: — ty^ for i — > TT, and it has the same behavior on A2 for bs — bA — 4, i.e. there 
occur scalar curvature singularities at the points C or D.) 

In order to obtain the metric potentials on T-L{ in terms of the Ernst potential and 
the constant Qq, we replace the potential u by a potential v in a neighborhood of l-Lf via 

Then we find 

nr. Q = -, ^ aim, e^('^'^) = c^, (108) 

flo Jf 

where the integration constant c in the expression for M can be determined from 

continuous transition to the axes. 

In the present case 6^ 7^ 6^, it turns out that "Hf is a regular Cauchy horizon, 

generated by the linear combination dp^ — aodp^ of the two Killing vectors. Regular 

coordinates can be introduced via 

TT — t — arcsin Vf, pi — p[ + ^ In t' p2 = P2 + In t', (109) 

where ki and /t2 are two constants that can always be chosen such that the resulting 
metric potentials are regular. 

Finally, we may look again at the singular cases 6^ = 6^ ±4. As mentioned earlier, 
the corresponding Ernst potential and the Kretschmann scalar on Ai (for 6^ = 6^ + 4) 
or A2 (for bs — bA — 4) diverge in the limit i — > tt. Since we therefore cannot find a 
solution of the LP on T-Lf that is continuously connected to the axes, it is not possible 
to construct the Ernst potential on Hi in these two singular cases directly. However, 
in order to study the situation on H{ in these cases too, we can consider a sequence 
of solutions with 6^ 7^ 6^ ± 4. Then, for each element of the sequence, the LP can be 
solved along all four boundaries of the Gowdy square and the corresponding expression 
for the Ernst potential S{ on Jif, constructed from this solution, is valid. It turns out 
that the limit i — >■ tt of £f remains regular for < ^ < tt, whereas Sf diverges as expected 
at C or D. Hence we can conclude that only the boundary points C or D of Hf become 
singular and the interior of T-Lf is still a regular null hypersurface. 



5. Discussion 



In this paper we have studied the class of smooth Gowdy symmetric generalized Taub- 
NUT solutions as interesting examples of Gowdy spacetimes with spatial topology. 
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This class is characterized by a special behavior of the metric potentials in a vicinity 
of the initial surface Hp {t = 0) which, in particular, implies that Hp is a smooth 
(past) Cauchy horizon. Utilizing Fuchsian methods, we were able to show that for 
smooth asymptotic data, describing the spacetime at Tip, there always exists a unique 
smooth Gowdy symmetric generalized Taub-NUT spacetime as solution to the Einstein 
equations for t e (0,7r). In a second step, we have investigated the behavior of these 
solutions on the symmetry axes Ai {6 = 0) and A2 {9 = -k). Using the complex Ernst 
formulation of the field equations and its reformulation in terms of an equivalent linear 
problem, we have constructed explicit formulae for the metric potentials on Ai and 
A2 in terms of the data on "Hp. Afterwards, it was possible to extend the solution to 
the future boundary {t = tt) of the Gowdy square and to find explicit expressions 
for the metric potentials there, too. It followed from these expressions, that we have to 
distinguish between four types of asymptotic data, which are characterized by the values 
Ba and &b of the imaginary part b of the Ernst potential at the points A {t = 6 = 0) 
and B {t = 0, 9 = 7i) and which lead to solutions with a completely different behavior 
on T-Li: 

(i) bB = bA + 4: 

In this scalar curvature singularity occurs at point C (t = n, 9 = 0). 

(ii) bB = bA- 4: 

Here, a scalar curvature singularity occurs at point D {t = 9 = n) . 

(iii) bs = bA- 

The spacetime is regular in the entire Gowdy square. In particular, the Ernst 
potential S and the metric potentials u, Q and M are smooth. Moreover, "Hf is a 
smooth Cauchy horizon, generated by the Killing vector (9p^. 

(iv) bB ^ bA and bs ^ bs ± 4: 

The spacetime is regular in the entire Gowdy square and the Ernst potential S is 
smooth, however the metric potential u is not well adapted to describe this case 
and blows up at Tif. However, there is no physical singularity at Tif. Instead, Tif 
is a smooth Cauchy horizon, generated by the null vector dp-^ — a^dp^ . 

This shows that — with exception of the two singular cases (i) and (ii) — smooth 
Gowdy symmetric generalized Taub-NUT solutions (with a past Cauchy horizon at 
t = 0) always develop a second Cauchy horizon at t = vr. This future Cauchy horizon, 
in the same way as the one in the past, is homeomorphic to and its null generator has 
closed integral curves. Hence our results can in particular be understood as a partial 
resolution of a problem which remained open in [10]. Namely, at least in our class of 
spacetimes, the union of the Gowdy square and its interior is isomorphic to the closure 
of the MGHD of corresponding Cauchy data. 

It is interesting to compare these results with the situation of spatial §^ x §^ topology 
as investigated in [19]. For this case, it was shown that Gowdy spacetimes with a regular 
past Cauchy horizon Hp develop a regular future horizon Hf if and only if a particular 
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quantity J, which can be read-off from the asymptotic datallj, does not vanish. In the 
hmit J — 7- 0, Tip transforms into a scalar curvature singularity. Hence, the behavior is 
similar to the S'^ case: with exception of singular cases, spacetimes with a past Cauchy 
horizon generically develop a future Cauchy horizon. However, the nature of the singular 
cases is slightly different: In the x §^ case, the curvature blows up along the entire 
future boundary Hf, whereas we find only singularities at the isolated points C or D on 
T-Lf for S'^ symmetry. 

Do our assumptions rule out important cases? This is not clear and probably 
difficult to answer. For example our assumptions do not allow solutions with past 
Cauchy horizons ruled by non-closed generators. There is no reason why such solutions 
should not exist. Another possibility excluded by our assumptions are solutions with a 
non-compact or incomplete Cauchy horizon. Indeed such solutions can be constructed, 
in the polarized case, by the techniques of Moncrief and Isenberg in (20] . 




Can our results be generalized to situations with less symmetry and eventually 
maybe even to generic solutions with Cauchy horizons? We have employed two, in 
principle, independent techniques for the two main steps of our discussion: the Fuchsian 
method for the basic existence proof and the soliton method for the study of the global 
properties of the solutions. As far as the Fuchsian method and the underlying singular 
initial value problem and hence the existence and uniqueness of solutions with prescribed 
"data" on a Cauchy horizon is concerned, we can say the following. In general, it 
cannot be expected that such an "initial value problem" for equations of hyperbolic 
type is well-posed. It seems at least necessary that the generator of the horizon, being 
a null hypersurface, is a Killing field. This is the case here and also in the more 
general f/(l)-symmetric case discussed by Moncrief [21]. Therefore under Moncrief 's 
assumptions, our Fuchsian technique should apply in the smooth f/(l)-symmetric case. 
The results in [231 [13] suggest that there must be a t/(l)-symmetry in a neighborhood 
of a Cauchy horizon in general vacuum spacetimes, at least if the horizon is compact 
and the generator has closed integral curves. So, there is hope that a similar singular 
initial value problem can be formulated under quite general assumptions and that the 
existence and uniqueness proof based on Fuchsian methods goes through. We also want 
to mention that most of the results here can easily be generalized to solutions with only 
finitely many derivatives, cf. Proposition 13.21 

As far as the Ernst formulation and the associated linear problem are concerned 
there is probably little hope of a generalization to situations with fewer symmetries. 
The introduction of the complex Ernst potential relies essentially on the existence of 
two Killing vectors, cf. (15T|) . (l53i) . and the linear problem makes use of the special 
structure of the Ernst equation. However, it should be quite straightforward to apply 
the methods to Gowdy symmetric spacetimes with additional electromagnetic fields. 
In that case one has to study the coupled system of the Einstein-Maxwell equations. 

In terms of the Ernst potential at Hp, J is defined as J = —-^^{bA ~ bs — 4Qp), where Qp denotes 
the constant value of the metric potential Q on Hp in the S'^ x §^ case. 
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for which, remarkably, a complex Ernst formulation and an associated linear problem 
exist as well. The corresponding calculations would follow closely the investigation 
of axisymmetric and stationary black hole spacetimes with electromagnetic fields as 
presented in pTSj . 
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Appendix A. Theory of second-order hyperbolic Fuchsian equations 

Here is a quick summary of the theory of second-order hyperbolic Fuchsian equations 
outlined in [6]; a more detailed presentation can be found in [7]. 

A second-order hyperbolic Fuchsian system is a set of partial differential equations 
of the form 

D^v + 2ADv + Bv -t^K^Av = f[v], 

in which the function f : (0, 5] x f/ — )■ M" is the main unknown (defined for some 6 > 
and some spatial domain U), while the coefficients A = A{x), B = B{x), K = K{t,x) 
are diagonal n x n matrix-valued maps and are smooth in x G f/ and t in the half-open 
interval (0, 6], and / = f[v]{t, x) is an n- vector- valued map of the following form 

f[v](t,x) := f(t,x,v(t,x),Dv(t,x),tK(t,x)d^v(t,x)y 

Here, we can assume that K is the identity matrix and that U is S^; the condition 
that [/ is a one-dimensional periodicity domain in the references above can easily be 
generalized. Of particular convenience for this is that the equations are geometric wave 
equations with spatial domain We assume that the time variable t satisfies t > 
and use the operator D := tdt to write the equations. We denote the eigenvalues of 
A and B by a^^\ . . . , a*^"^ and h'^^\ . . . , respectively. When it is not necessary to 
specify the superscripts, we just write a, h to denote any eigenvalues of A,B. With this 
convention, we introduce: 

Ai := a + \/a^ X2 := a — y/a^ — h. 

It turns out that these coefficients, which might be complex in general, are important 
to describe the expected behavior at t = of general solutions. 

Consider a second-order hyperbolic Fuchsian system with coefficients a, 6, Ai, A2. To 
simplify the presentation, we restrict attention to scalar equations (n = 1) and shortly 
comment on the general case in the course of the discussion. Fix some integers l,m > 
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and constants a, 5 > 0; indeed a can also be a smooth positive function on U. For 
w G C''{{0, 5], H"^(U)), we define the norm 



iY.J2 [ t^^"^^'^"^""^ \d',DPw{t,x)f dx 



0<t<<5 . „ „ 
— \p=0 9=0 

and denote by Xs^a,i,m the space of all such functions with finite norm ||w||5^Q,^/^m < oo. 
Throughout, H^{U) denotes the standard Sobolev space. To cover a system of n > 1 
second-order Fuchsian equations, the norm above is defined by summing over all vector 
components with different exponents used for different components. Recall that each 
equation in the system will have, in general, a different root function A2. We allow that 
a = {a^^\ . . . is a vector of different positive constants for each equation. The 

constant 6, however, is assumed to be common for all equations in the system. The 
motivation for including the quantity A2 into the definition of the norms is the standard 
singular initial value problem in [6]. Throughout it is assumed that 3fJA2 is continuous 
and it is then easy to check that {Xs,a,i,m, \\ ■ \\5,a,i,m) is a Banach space. For each non- 
negative integer / and real numbers 5,a > 0, we define X^^,; := C\p=o-^s,a,p,i-p, and 
introduce the norm 

/ ' X 1/2 

||/||5,Q,« := f X] 11/11 |a,p,i-pj , / e Xs,a,l- 

It turns out that we must also use spaces [X^^a^i-, || ■ IliTai)- These are defined as before, 
but in the norm of some function /, the highest spatial derivative term 9^/ 

is weighted with the additional factor t (in our case here). It is easy to see under 
the earlier conditions that also [X^^a.i-, || ■ ||7a«) Banach spaces. Let us also define 

Let US choose a function u on (0, 5] x f/ (whose regularity is fixed in the main 
existence theorem). The singular initial value problem associated with u is then defined 
as follows. We ask whether there exists a solution v of the given second-order hyperbolic 
Fuchsian system so that the remainder 

w{t, x) := f (t, x) — u(t, x), 

can be interpreted as "higher order" in t at t = 0, where u is interpreted as the leading- 
order term. By this we mean that w is an element in Xs^a,i for some (sufficiently 
large) a > 0. Often u will be parametrized by certain free functions which we call 
asymptotic data. An example of a singular initial value problem with a leading-order 
term parametrized by asymptotic data is Proposition 13.21 In that case Ai is zero for 
both equations and A2 is zero for the first equation and —4 for the second one. 

The existence and uniqueness result of the singular initial value problem is 
Proposition 3.5 of [B]. The existence proof is based on a new approximation scheme 
which can also be implemented numerically. See in particular [8| for more details. 
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